Impact of long-range interactions on the disordered vortex lattice 
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The interaction between the vortex lines in a type-II superconductor is mediated by currents. In 
the absence of transverse screening this interaction is long-ranged, stiffening up the vortex lattice 
as expressed by the dispersive elastic moduli. The effect of disorder is strongly reduced, resulting 
in a mean-squared displacement correlator (it 2 (R, L)) = ([u(R, L) — u(0,0)] 2 ) characterized by a 
mere logarithmic growth with distance. Finite screening cuts the interaction on the scale of the 
London penetration depth A and limits the above behavior to distances R < A. Using a functional 
renormalization group (RG) approach, we derive the flow equation for the disorder correlation func- 
tion and calculate the disorder-averaged mean-squared relative displacement (tt 2 (R)) oc ln 2<T (i?/ao). 
The logarithmic growth (2a = 1) in the perturbative regime at small distances [A.I. Larkin and 
Yu.N. Ovchinnikov, J. Low Temp. Phys. 34, 409 (1979)] crosses over to a sub-logarithmic growth 
with 2a — 0.348 at large distances. 



I. INTRODUCTION 

Disordered elastic systems have attracted much atten- 
tion in the last decade. They have applications in vari- 
ous physical systems, such as charge density waves j» do- 
main walls in Ising-ferromagnets^and vortices in type-II 
superconductors^ The most interesting physical proper- 
ties of such disordered elastic systems are their struc- 
tural (static) order and their dynamical response under 
external forces. The dynamics involves a finite criti- 



cal force with cree 



4,5,6 



at small- and depinning"* at 



large forces. Here, we concentrate on the structural as- 
pects as expressed through the displacement correlator 
(u 2 (r)) = ([u(r) — u(0)] 2 ) (angular brackets denote av- 
erage over disorder). In general, the mean-squared dis- 
placement grows with distance as (u 2 (r)) oc |r| 2 ^, where 
£ is the wandering exponent characterizing the roughness 
of the elastic manifold; for £ > the disorder is relevant 
and the manifold is roughs 

The generic model describing the competition be- 
tween elasticity, quenched disorder, and possibly ther- 
mal fluctuations involves an isotropic homogeneous elas- 
tic medium (manifold) characterized by d internal dimen- 
sions (r) and an iV-component vector field for the dis- 
placement u(r) ([d + N] -model). The system we consider 
in the following is the flux-line lattice in a type-II super- 
conductor submitted to quenched disorder. The distor- 
tions of the lattice are described by three internal dimen- 
sions r = (R, L) and a two-component displacement field 
u(r) ([3 + 2]-model). 

A closer inspection shows that the underlying physical 
properties of the vortex lattice infer interesting modifica- 
tions of the generic structural characteristics, both quan- 
titatively and even qualitatively (a similar type of spe- 
cific properties is found in other systems, e.g., in Wigner 
crystals) . The underlying structure of the vortex lattice 
introduces a number of additional length scales in the 
problem: i) On a microscopic level, the flux-line lat- 
tice is built from interacting lines; their arrangement 
in a periodic fashion with distance ciq drastically re- 



duces the effect of disorder on large scales a) The 
line-nature leads to a tensorial elasticity with the shear 
modulus different from the compression and tilt mod- 
uli; this introduces a weak dependence of the wander- 
ing exponent upon the ratio of compression and shear 
moduli^ Hi) Long-range interactions between the vor- 
tex lines stiffen up the lattice as expressed by dispersive 
compression and tilt moduli^i^ thereby qualitatively 
reducing the effect of disorder; in the presence of trans- 
verse screening, this effect is limited to scales below the 
London penetration length A. iv) The long-range interac- 
tion between the vortices and the defects implies a long- 
range correlation in the disorder landscape, modifying 
the physics of thermal depinningAi& 

In this paper, we focus on the effect of the long-ranged 
interaction between the vortices, see Hi) above. This 
interaction plays an important role in the short- to in- 
termediate distance regime before being cut off by trans- 
verse screening at R ~ A in transverse-, L ~ A 2 /ao in 
longitudinal direction. Let us put the results for the dis- 
placement correlator into context and follow its evolution 
from small to large scales: 

Because of the line-nature of the vortex lattice (cf. ii) 
all results are anisotropic; we will discuss the limits in 
the plane (R) and along the direction of the vortex lines 
(L) separately and comment on the general behavior as 
appropriate. We assume weak point-like disorder; the 
effect of the interaction changes the algebraic growth of 
the displacement correlator (it 2 (r)) into a mere logarith- 
mic growth, as can be seen in lowest-order perturbation 
theory^ valid at small distances R < R c and in the weak 
pinning limit , 



In [R/aao] 
In [R c /aa ] 



(1) 



Here, £ is the correlation length of the disorder land- 
scape properly defined in Eq. I|41|l (for point-like disorder 
the scale £ is of the order of the vortex core dimension) 
and a is a number of order unity. Later, we will obtain 
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a w c/y/w, where c = e~ 7 and 7 w 0.577 is the Euler con- 
stant. The result (u 2 (L)) in the longitudinal direction is 
obtained from the above equation via the substitution 
R — > [(c/y / 7r)aoi] 1 ^ 2 . The dependence on the disorder 
strength has been encoded in the collective bundle pin- 
ning radius R c defined via (u 2 (R c )) — £ 2 ; for suitably 
weak disorder we have ao < R c < A and we will discuss 
other cases later. The long-range nature of the interac- 
tion between vortices is relevant for 



ao < R < A, ao < L < A 2 /ao- 



(2) 



When (u 2 (r)) increases beyond £ at R c the appearance 
of multiple minima leads to a breakdown of the perturba- 
tive analysis For R c < R < A, we use the renormaliza- 
tion group (RG) to calculate the mean-squared displace- 
ment and find the result 



(u 2 (R))~e 



hx[R/a ] \- 
ln[R c /a }J 



(3) 



exhibiting a sub- logarithmic growth with a « 0.174 (to 
the precision assumed here, we drop the numerical a from 
Eqs. iJSJ to J5J); the substitution R — > y/a L again pro- 
vides the result along the longitudinal direction. Equa- 
tion J3J is the central result of this paper and will be 
derived below. 

On larger length scales, the static properties of the 
disordered flux-line lattice are well understoodAi^ At a 
distance R a the displacements become of the order of 
the lattice constant a , (u 2 (R a )) ~ a 2 , and the effect of 
disorder is reduced since the flux-lines no longer probe in- 
dependent disorder realizations. In the absence of trans- 
verse screening the correlator then has an even slower 
growth (u 2 (R)) oc lnln[i?/i? a ]. 19 In a real system, how- 
ever, the vortex-vortex interaction is usually screened on 
scales A smaller than R a , the lattice gets softer and there- 
fore rougher. For A < R < R a (random manifold regime) 
one obtains^ 



(u 2 (R))~£ 



( g/A 



Vln[i? c /a ] 



2C 



(4) 



The result in longitudinal direction can be obtained via 
the substitution R — > clqL/X (the anisotropy is modi- 
fled as compared to the case R < A above). The func- 
tional RG calculation^ provides a value of the rough- 
ness exponent £ which depends weakly on the ratio of 
the compression (en) and shear (cqq) moduli, assuming 
the value £ ~ 0.174 in the limit en 3> cqq relevant in the 
vortex lattice. Finally, periodicity becomes important at 
R a (Bragg Glass regime) and one recovers a logarithmic 



(u 2 (R)) 



In 



R_ 

Ra 



(5) 



where the prefactor of the logarithm does not depend 
on the disorder strength (the latter is encoded in R a ). 



The physical properties of the vortex lattice, its period- 
icity and the long-ranged interaction between the con- 
stituents, thus qualitatively reduce the effect of disorder. 

The paper is organized as follows: In the next section 
we present the model and perform the perturbative cal- 
culation which provides the starting point for the renor- 
malization group analysis. Section III is devoted to the 
derivation of the flow equation for the disorder correlator 
and the calculation of the effective mean-squared relative 
displacement. We summarize and conclude in section IV. 



II. THE MODEL 

We consider an isotropic type-II superconductor in the 
mixed state and use a coordinate system with the mag- 
netic field B applied along the z-axis. In the absence of 
disorder the vortices form a triangular Abrikosov lattice 
with a lattice constant a A , a\ — 2a 2 ,/v3, and a 2 , = &o/B 
($0 = hc/2e is the magnetic flux quantum). The disor- 
dered vortex system is described by a free energy involv- 
ing both clastic and pinning parts. Within the continuum 
elastic theory, the elastic part takes the forrn^ 



T c \ ful = - 



d 3 k 



[ Cll (k)|K-u(k)| 2 



2J (2tt) 3 

+c 66 |K ± •u(k)| 2 + c 44 (k)fc 2 |u(k)| 2 ] 



(6) 



where k = (K,k z ) = (k x ,k y ,k z ), K± = (k y ,-k x ), and 
with the integration defined over the lattice Brillouin 
zone. The compression and tilt moduli are strongly dis- 
persive, i.e., they depend on the wave vector k of the 
distortion, 



cn(k) w c 44 (k) 



c 44 



1 + A 2 fc 2 



C44 



B 2 
4tt' 



(7) 



where k = |k|. The dispersion is due to the long-range 
current-mediated interaction between the vortices and is 
relevant in the intermediate distance regime, 



ao < R < A, a < L < A 2 /ao 



(8) 



(r = (R, L), R = |R|), before it is cut on the scale of 
the penetration depth A; we will restrict ourselves to this 
intermediate regime in the following. As usual, we ignore 
the weak dispersion in the shear modulus c 66 , 



C66 



(8ttA) 2 



(9) 



The static Green function corresponding to the elastic 
energy JBJ reads (/i, v = x,y) 



G^(k) 



*T(K) 



Cii(k)X 2 + C4 4 (k)fc2 



(10) 



C 6e K 2 + C44(k)fcf 
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with the longitudinal and transverse projection operators 
V^(K) = K^Ky/K 2 and V^(K) = 8^ - V^{K), re- 
spectively (K = |K|). In the regime we are interested in, 
the compression modes involve large energies (en 3> Cqq) 
and the longitudinal part of the Green function can be 
neglected; we denote the transverse part by G^ . 

Assuming weak point-like defects, we describe the dis- 
order term in the free energy via the pinning energy den- 
sity -Ep in (r,u(r)), 



jF pin [u] = J d 3 r £ pin (r,u(r)). 



(11) 



The pinning energy density -Epi n derives from the con- 
volution of the pinning potential with the vortex form 
function 3 and is characterized by the disorder correlator, 

(£p in (r, u)£ pin (r', u')> = S 3 (r - r')K (u - u'), (12) 

where K is proportional to the disorder strength and is a 
short-range function decaying over the characteristic size 
£ of the vortex core. For the disordered vortex lattice the 
correlator has an algebraic decay ^> Kq cx (1/u 2 ) In (it/£) 
for u < A. We consider weak disorder where the vortices 
are pinned in bundles of size R c > ao; this contrasts with 
single vortex pinning realized for larger disorder strength 
to be discussed later. The above description remains 
valid as long as the disorder-induced displacements re- 
main small, (u 2 ) < a 2 ,, such that the vortices probe a 
different random environment upon deformation. As the 
displacement becomes of the order of a on large scales, 
one has to take the periodicity of the lattice into account. 

The equilibrium configuration u(r) of the system is 
obtained by minimizing the free energy functional .F[u] = 



pin and is implicitly given by 



u"(r) 



d 3 r'G^(r-r')^E^(r',n^% (13) 



where (r) is the Fourier transform of the Green func- 
tion I|1U|) . This 'equation of state' allows for calculating 
the mean-squared displacement correlator (u 2 (r)). 



A. Perturbation theory 

For small displacements (u 2 (r)) < £ 2 , the u depen- 
dence in the force (— d u «E pm ) can be neglected in (jT3|l . 
and the mean-squared relative displacement can be writ- 
ten in the form 



(u 2 {v)) 



d 3 k 

xG w (k)G pi/ (-k), 



(1 — coskr) 



(14) 



with the integration restricted to the lattice Brillouin 
zone in the plane. 

We denote derivatives of the correlator with respect to 
u v by a superscript, 



and sum over indices that appear twice. Equation 114|) 
simplifies since K$ is isotropic: Kq x (0) = Kq V (0) = 
Kq(0) and Kq V (Q) = (the primes denote derivatives 
of Kq(u) with respect to u = |u|). Inserting the trans- 
verse part of the Green function (|10f) . we obtain 



(u 2 (r)) » -2X '(0) 



d 3 k 



1 — cos kr 



(2tt) 3 [c 66 K 2 + c M (k)k 2 ] 2 - 



(15) 



In the intermediate regime, ao < R < A, ao < L < 
A 2 /ao, the relevant values of k z are much smaller than 
K. Hence, we can approximate c 44 (k) = c 44 /(l+A 2 fc 2 ) w 
C4,i/\ 2 K 2 and obtain 



(u 2 (v)) « -2^(0) 



d 3 k [1 - cos(kr)]A 4 J ft: 4 
(2tt) 3 {c e6 \ 2 K 4 + c 44 k 2 z }< 



(16) 



The divergence of the Green function as k — * is cut 
off by the numerator on the scale k ~ 1/r. For large 
wave vectors, the integral is limited by k ~ 1/clq (bundle 
pinning regime). In a first analysis, we therefore replace 
CEIi by 



(u 2 (v)) 



^o'(O) 1^ J d 2 Kdq- { 



K 4 



[K 4 + q 2 

where q = k z y/ 'c 44 / 'cqq / 'A « /c z 4y / 7r/ao and 



(17) 



(18) 



47TC 66 V / C66C44 

The integral in (|17|l is essentially four-dimensional, 

2 K 4 d 2 Kdq ^ 1 d 3 (K 2 )dq 1 d 4 y ^dy 

^ [K 4 + q 2 } 2 [(K 2 ) 2 + q 2 } 2 ~ 2vr 2 y 4 ~~y~ 

(the symbol = reminds that we are using d d x — 
x d ~ x Sddx, with Sd the surface of the d-dimensional 
sphere), and has a logarithmic divergence at small y = 
(K 4 + q 2 ) 1 / 2 . The integral l|17|l is dominated by values 
q r~j K 2 ; for R 2 ^> clqL the small scale cut-off is given 
by y ~ 1/-R 2 , while a cut-off y ~ 1/clqL applies in the 
opposite limit L 3> R 2 /qq. We therefore obtain 3 - 



( U 2 (r))«-K '(0)/ln 



P(r) 
ao 



(19) 



where the function p(r) assumes the limits 

p(R 2 > a L) = R/a , 
p(R 2 ^a L) = (a Q L) 1/2 /a', 



(20) 



with a and a' numbers of order unity; the proper inter- 
polation for arbitrary values of R and L is beyond the 
accuracy of the present analysis. Note that as a conse- 
quence of the dispersion in the tilt modulus, the lengths 
along the longitudinal (L) and transverse (R) directions 
scale differently, L ~ R\J c^j c m ~ R 2 /a for R < A, 
L < A 2 /ao. The logarithmic growth of (|19fl is a conse- 
quence of the dispersion in the tilt modulus (and hence of 



4 



the long-range interaction) , effectively lifting the problem 
to four dimensions where the disorder is only marginally 
relevant the prefactor of the logarithm is proportional 
to the disorder strength. 

The simplifications made in 117fl do not allow for a de- 
termination of the numericals a, a', and the shape /?(r); 
in the following we perform a more accurate calculation 
of the original expression i|14|) in order to obtain a pre- 
cise result for the argument under the logarithm in (|19|l . 
We approximate the hexagonal lattice Brillouin zone by 
the circular version K < K BZ = ^/in/ao; using cylindri- 
cal coordinates k = (K, 0, fc z ), the angular integration in 
(|16|l can be carried out directly and the integration over 
the fc z -axis can be performed using residues, after which 



2 (r)> 



Abz 



IKg(0) 



1 - T n (KR)e~ a ° LK / 4 v^ 

d.K—^% . (21) 

K 



For R or L are much larger than do, the integral in 1211) 
can be approximated by 



A{R, L) 



.1 



-a LK 2 



dK- Jo{KR)e ~ e-Wez, 
K 



where the constant c is fixed by demanding 
" Kbz dK f°° dK 



K 



K 



cK'/Ki z 



0.56, where 7 w 0.577 



for (5^0; one obtains c = e~ 
is the Euler constant. 

The integral A(R,L) exhibits a logarithmic behavior 
for small K which then is cut off cither by R or by 
\/Lciq. The approximation of a circular Brillouin zone 
(see above) is correct up to a factor of order unity under 
the logarithm. 

We directly evaluate the integral A(R,L) and obtain 



A(R,L) = In 



R 



(c/0r)ao 



1 



r 0, 



irR 2 



(r(0, x) — J^ a dte^ t /t is the incomplete gamma func- 
tion), displaying a complicated dependence on the cut- 
off parameters R and L. In the limit of large R, R 2 3> 
2ca 2 /ir + 2aoL/y/ir, the above result simplifies to 



A(R,L) w In 



R 



(22) 



while in the opposite limit of large L, R 2 <C 2ca\j-K 
2qqL/ 'yV (and L 3> caof-y/n) one obtains 



A(R, L) « In 



L 



(c/V7r)a 



1/2 



(23) 



The present, more accurate, analysis then allows for the 
definition of the cut-off function p(R, L) for arbitrary val- 
ues of R and L within the perturbative regime, 



Comparing the above results with (|19(l and (|20ll , we find 
the numericals a « c/y/n w 0.32 and a' = a 1 ' 2 « 0.56. 
Defining the collective pinning radius R c via 



(25) 



we obtain Q with the given value of the numerical a. 

The above perturbative results are valid within the col- 
lective pinning regime R < i? c , which is equivalent to 
p(R,L) < Rc/a ~ 3.16 R c . Beyond this regime, the per- 
turbative result l|19[) breaks down because of the appear- 
ance of multiple minima as the displacement (u 2 (r)) 1 / 2 
increases beyond the characteristic scale £ of the disorder 
potential. 17 Along the longitudinal direction, this break- 
down occurs at the larger scale Lj? 3A6R 2 /ao > R c - 
Within the present perturbative approach the length 
scale £ appearing in 125|) has been introduced ad hoc — 
this will be different in the functional RG treatment ex- 
tending the discussion to larger length scales beyond R c . 



III. FUNCTIONAL RG EQUATION 

The logarithmic behavior of the displacement correla- 
tor (|19f) provides the motivation for applying the renor- 
malization group 20 (RG); we briefly summarize the main 
steps in the derivation of the functional RG equation for 
the disorder correlator Kq(u) using the real-space renor- 
malization procedure introduced in Ref. l2lL We aim at 
finding the effective disorder correlator on large scales 
(R c < R < A, £j? < L < A 2 /a ). The elastic coefficients 
are not renormalized as the pinning part of the system 
is invariant under the 'tilt' u(k) — > u(k) + v(k) for any 
vector-function vm^^ 

Assuming renormalization up to a distance (Rx,L\), 
we define the displacement ui, the energy E^ n , and the 
correlator K\ at the scale (i?i, Li). We perform the RG 
step to i?2 > Ri, L2 > Li, and split Ui = U2 + w (using 
the self consistency equation Q13[)) into a far- 

<(*) = -[ dVG^(r-0^i5$(iWr / )) 
Jn> 

with Q> = {|R - R'| > R 2 ; L - L' > L 2 }, and a near- 
field contribution, 

«/» = - / d 3 r' G^(r-r')d u ,E^l(v', Ul (r')), 
Jn 

where Q = < |R- R'| < R 2 ; L t < L - V < L 2 }. 
We express the effective free energy E^(u 2 (r), r) for the 
far-field contribution as a Taylor expansion in the near- 
field w(r) and, assuming Gaussian disorder, determine 
the effective pinning energy correlator 

K 2 (u 2 - u' 2 ),5 3 (r - r') = (^(r, u 2 ) E^r', u 2 )) 



p(R,L) 



a e 



A(R,L) 



(24) on the scale (R 2 ,L 2 ). Evaluating the terms to second 
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order in X\ (first loop), one obtains 

1 



K 2 (u 2 )=if 1 (u 2 )+/^ 
where 



(26) 



j^ pK = A d 3 r / qh"(j. _ r ')GP K (r - r'). (27) 

The Green function is most conveniently expressed in 
Fourier space, 



d 3 k 



G^(k)G pK (k), 



(28) 



where J7' is the appropriate &-space integration domain 
corresponding to fl (see below). Since en > cq§ we can 
neglect the longitudinal part of the Green function I|1Q|I 
and inserting the transverse part Gt into (|27|1 the ten- 
sorial structure of the integral can directly be calculated 
using cylindrical coordinates (X, 4>,k z ), 

/ d<j)T^(K)T^(K) = — A^ pK , (29) 
Jo 8 



where 



(30) 



(for the isotropic situation with G^ v oc 8^ one obtains 

Equation (|28|l resembles the equation for the displace- 
ment correlator l|14|) ; following Ref. [2J we approximate 
the term 1 + X 2 k 2 in c 44 (k) by X 2 X 2 and use the fact, 
that the resulting integral is effectively four-dimensional, 
cf. (|17|l . Performing the RG-step only in transverse direc- 
tion (L 2 — Li), the /^-integration extends over the whole 
axis, while the if- integration is cut-off at if 2 ~ 1/-R 2 and 
Ki ~ l/i?i; in the opposite case, the integral is cut by 
k Z 2 ~ 1/^2 and k z ± ~ 1/L\. The final result is conve- 
niently expressed^ through the function /?(r) as defined 
in |2U|) 



16^ m vp(^i,ii: 



(31) 



note that, given the accuracy of the present analysis, we 
do not make use of the result (|24|l which has been derived 
within the perturbative regime only. We define the flow 
parameter 



l = Hp(R,L)/a ], 



(32) 



where the starting point of the flow is chosen by matching 
our RG analysis with the result from perturbation theory, 
cf. PI). 

Inserting (|3I|I into (|26[) yields the RG equation for the 
effective pinning energy correlator, 



f 



diKt(u) = —A^ K ( -K^Kr - K? K (0)K 



*up 



(33) 



Summing over pairs of indices, this simplifies to 



I 



diKi = [ K^Xr - 2X^X^(0) + 



(34) 



The first line of this RG equation contains the terms 
that are obtained for an isotropic situation2&2£ with 
G pv oc 8^ v ; the terms in the second line arise due to 
the transverse structure of the Green function entering 
l|29p . Up to scaling terms, an equation equal to (|34|l is 
obtained in the non-dispersive [3+2]-case for en 3> Cqq, 
while a different tensorial structure is found assuming 
arbitrary values of en and cqq, cf. Ref. Il3l 

In the regime we are interested in, the correlator is 
isotropic, Xi(u) = Ki(u), and the derivatives can be 
rewritten as 



Xr(u) - 
X^(0) = X[\Q)5^- 
the RG equation l|34|l then takes the form 



(35) 



diKAu) = — 
16 



-4A',"(o)(k;' + S 

V u 



(36) 



where the primes denote derivatives with respect to 
u = |u|. Taking four derivatives with respect to u and 
assuming an analytic correlator, one obtains 

fl^(0)=J[xf 4 >(0)] a . 

At the scale l c , the fourth derivative of the correlator 
becomes infinite at the origin, Kjc^ (0) = oo, and the 
second derivative develops a cusp singularity. This non- 
analyticity provides a definition for the collective pinning 
radius R c which depends on the fourth derivative ifg 4) (0) 
of the initial correlator, 



lc. = 



1 



ii^(0) 



In [Rc I aa ] 



(37) 



Next, we determine the static structure of the vortex 
lattice on large scales. As for the disorder correlator Xi, 
we use perturbation theory to find the renormalized dis- 
placement field (u 2 (r)) on scale r. We split the integral 
(|17|l into shells and replace X^(0) in each shell by the 
scale-dependent disorder correlator K"(0). Integrating 
over all shells we obtain 



(u 2 (v)} 



Anp(r) 

I / dl'K'/,(0) 

J In an 



(38) 



As in (|17J) . the integral is cut off by the function p(r), 
which is defined in (|2(Jf) via its limits for R 2 much larger 
or much smaller than a^L. 
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IV. ANALYSIS OF THE RG FLOW 

Analyzing (|36l) . the disorder correlator K\ flows to- 
wards zero on large scales, indicating that the disorder 
is only marginally relevant. Here, we are interested in 
learning how K','(0) flows to zero with increasing I in or- 
der to integrate the equation (|38l) for the displacement 
correlator. We will see that a proper rescaling of the 
displacement field u(r) and of the correlator Ki maps 
the flow equation (|36|l describing the marginal situation 
to the flow equation in one dimension Ies9i2*2&; this will 
allow us to make use of results derived from an e = 4 — d- 
expansion at the value e = 1. 

Following (|38[1 we need to integrate the flow equation 
for K'/iO). Differentiating (|34|l or l|36l) twice with respect 
to u and evaluating the result at u = one obtains 



(39) 



In the Larkin regime, R < R C1 L < Lj?, the disorder cor- 
relator is an analytic function (K' l "(0 + ) = 0) and hence 
K'i'<i (0) = (0) remains constant; integrating we 
confirm the perturbative result (|19J) . 



rKlc 

l (r))=I 

Jo 



dl'[-K' l ',(Q)] = -K^0)Il, (40) 



with I — ln[p(i?, L)/aao]. This result is in agreement 
with the statement in Ref. [23] that higher order terms 
in perturbation theory do not contribute to (u 2 (r)). 

Combining 1)3 7|) and 14U|) at p = p{R c ), we find the 
proper definition of the relevant length scale £ of the dis- 
order landscape in terms of derivatives of Kq, 



K, 



(4), 



£ 2 ; 



(41) 



this precise definition of £ naturally appears within the 
renormalization group framework. 

Going beyond l Cl it is more convenient to investigate 
the asymptotic flow of Ki(u), rather than solving Eq. 
(|39|l . Since the magnitude Ki(0) decreases and the func- 
tion Ki(u) broadens, we perform the rescaling 



16"'" l- Jt \l-J (42) 
in order to extract the characteristic geometry of the cor- 
relator (its height and width). The simplified flow equa- 
tion for /; then will converge to a regular and finite func- 
tion /* for large I such that we arrive at a proper asymp- 
totic scaling form for K\ ~ l~ T f*(u/l a ). Inserting lH^ 
into i)36|) , the resulting flow equation involves terms pro- 
portional to l~ T ~ 1 and Z~ 2r ~ 4<T ; we require these terms 
to scale the same way in order to achieve a flow /; — * /* 
and hence we fix r = 1 — 4cr. The flow equation for /; 
then reads 



(43) 



-/r4 + ^m _4//'(o) (//'+£ 



where the first two terms on the rhs arise due to the ex- 
plicit ^-scaling in (|42|l and the primes denote derivatives 
with respect to u = u/l a . As we will see, this equa- 
tion is the same (up to InZ — > I) as the (first loop) flow 
equation for the rescaled disorder correlator Ri in the 
non-dispersive [3+2]-case in an e = 4 — d-expansion, cf. 
Eq. (IZl below. 

In the dispersive [3+2]-case, the system is at the upper 
critical dimension and the relevant integral for I^ u p k t cf. 
(|27|l . exhibits a logarithmic behavior, cf. (|31|l . This leads 
to the flow equation (|33|l in the form d\Ki — 1 _ loop (Kf). 
In the non-dispersive case the upper critical dimension is 
fouri£ In d = 4 — e the integral (|27I31|) has an algebraic 
behavior 

^u P K = _L A ^ K - (p(v 2 y - p(n) e ) (44) 
16 e 

(with / = 21/ ir A for e = 1), leading to the flow equation 
dptKj = Oi_ loop (A'?), where the flow parameter is 



I = In 



A 



(45) 



(here, C44 rs 644 and therefore p{R, L) is different from 
(H3J): P(R,0) ~ R, P(0,L) ~ aoL/X). Next, we rescale 
the correlator Kj = exp(— el) Kj and obtain the flow 
equation in the form diKj{u) — eKj + C 1 _ loop (A'?). Fi- 
nally, we perform the rescaling 



ffiiu) 



e- al R r 



(46) 



corresponding to l|42l) above if we replace I — ► exp(Z); the 
latter accounts for the change from a logarithmic to an 
algebraic behavior when going away from the marginal 
dimension. Carrying out the above sequence of steps, 
we arrive at the standard flow equation^ (we define 
u = ucxp(— Q) and assume en ^> cqq; the primes de- 
note derivatives with respect to u) 



In, 
di 1 



(e-4C)Rj 



CuR'j - 



(47) 



R'r 3 (R'tY { R~ 

1 u 2 \ u 1 \ 1 u 



where we have used a — — 4£. For e = 1, this equation 
is identical to (|43|l if we replace dh\l — > dl and a — > £. 
Comparing to standard derivations of flow equations of 
the type ijlTI) . cf. Refs. |2^ and IT3L here, we have de- 
layed all the rescaling of the field (u) and of the coupling 
function (Kf) to the very end of the calculation. The 
rescaling (|46|) then again is motivated by the desire for 
a flow equation for Rt(u) admitting a fix-point function 
R*(u) corresponding to Q^ 2 w 0.174 as found in Ref. 
for en ^> C66 (we denote by Q L N the one- loop value of 
the wandering exponent for the [d + A]-model). 
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We briefly remind the arguments of Ref. [2j| provid- 
ing the proper fix-point function R*(u) and the wan- 
dering exponent £. Depending on the starting condi- 
tions, the numerical integration of the fix-point equation 
for R* provides three classes of decaying solutions: for 
C < Qri t ne function R*(u) crosses the M-axis at least 
once. Choosing £ > Q L 2 , the solution is always positive 
and decays in a power law fashion for large u. Finally, 
for £ = Q L 2 , R*(u) is positive as well and decays expo- 
nentially. The first class of solutions is forbidden by the 
flow as a positive initial correlator remains positive un- 
der the RG flow, hence £ > Q L 2 . On the other hand, any 
algebraic tail Rjiu) ~ Aj/ if vanishes under the flow, 



at large distances 



dl 1 



(1 - 4C - jOAj < (1 - 4^ 2 - iC^)A h 



provided that 7 > 1/Q L 2 — 4 (this is the usual power- 
law separating short- from long-range correlator a 23 ! 25 ); 
for our correlator with Kq cx 1/u 2 , we indeed obtain that 
the tail disappears under the flow and the only accept- 
able solution is the exponentially decaying short-range 
fix-point function realized for £ = Q L 2 w 0.174. 

Next, we can make use of the analogy between the 
dispersive four-dimensional and the non-dispersive three- 
dimensional case discussed above: the fix-point equations 
for R* and for /* are identical and we obtain an exponent 
a = Q L 2 « 0.174 describing the power-law type evolution 
of the correlator's width, cf. I)42[l. An important differ- 
ence appears when considering higher-loop corrections: 
such corrections to R* and the roughness exponent £ are 
relevant for the e-expansion, 28 ' 2 ^ since under the trans- 
formation (j46(l each higher-loop term appears with an 
additional factor exp[— (a + 4£) I] = 1. On the other 
hand, in the marginal situation the higher-loop correc- 
tions to // are irrelevant; higher-loop terms come with 
an additional factor ^( T + 4cr ) = ; _1 and thus can be ne- 
glected at large I. This happens because at the upper 
critical dimension the correlator Ki flows to zero, and 
higher order corrections in Ki flow to zero faster and can 
be neglected, whereas Kr (for e > 0) flows towards infin- 
ity, such that higher order terms may be important (the 
growth in amplitude and width turn out to compensate 
each other). The above one-loop value for a thus is exact. 
This is similar to what happens in critical phenomena, 
where the leading singularity can be calculated exactly 
in one loop at the critical dimension. 

We now are in a position to determine the displacement 
correlator (it 2 (r)) on large scales following l|38fl . Taking 
the second derivative of (|4^|l and using /; — > /*, we ob- 
tain K'^ lc (0) ~ (16//)/*"(0)H 1 - 2<T) . We assume that 
this asymptotic behavior crosses over to the solution at 
small distances at / ~ l c . Using K" (0) = Kq(0), we 
obtain the large I behavior 



^ ic (0)~^(0) 



1-2CT 



( U 2 (r)) ~ I / dl'[-K'/,(0)] 



~ -IK'^l ■ (48) 

Using l|41|) and l c — \n(R c /aao), we finally arrive at 

(u 2 (r)) ~ e(l4 4) (0)l) 2a (49) 

_ ^2 ( ln [p(r)/«o] 
\\n(R c /aa ) / 

with p given by l]20p. Such a sub- logarithmic behavior 
has already been conjectured in Ref. 0; here, we have 
put the sub-logarithmic growth on a firm basis and have 
accurately calculated the value of the exponent. 

The sub-logarithmic growth found in H49|) above has 
to be cut off as the interaction between the vortices is 
screened on distances beyond the penetration depth A or 
when the displacement field increases beyond the lattice 
constant a$. If the latter condition ((u 2 (i?)) 1 / 2 > ao) 
is realized within the dispersive regime the growth in 
the displacement field turns even slower: beyond u > ao 
the effective disorder correlator becomes a periodic func- 
tion in u, implying that a = in (|42|l and hence 
K" >x (0) oc 1/1; the integration in (|48l) then leads to a 
growth (u 2 (r)) cx lnln[/o(r)/ p(R a )]^. However, typically 
one would expect that screening becomes relevant before 
periodicity in Kq and the displacement correlator follows 
the evolution dictated by the usual non-dispersive ran- 
dom manifold- and Bragg Glass regimes, cf. (@J and JSJ. 

The above results have been derived within a weak col- 
lective pinning approach with R c > ao and are valid on 
scales p(r) > ao. On the other hand, increasing the disor- 
der strength, the collective pinning radius R c drops below 
ao and vortices become individually pinned. 3 Within this 
single vortex pinning regime, finite segments of length 
~ [^oAo'HO)] 1 ^ 3 are pinned by the collective action of 
defects; here, ko(u) denotes the pinning energy correla- 
tor for an individual vortex line and So = ($o/47rA) 2 is 
the short wave-length line tension. Given the same un- 
derlying disorder potential, the single vortex- and vortex 
lattice correlators are related via k (u) = aoK (u). The 
collective pinning length L c and the collective pinning 
radius R c then are formally related via 



Rc. 



ao exp 



const. I — 
a a 



(50) 



Integrating (|38|l . we obtain the displacement correlator 



cf. H37fl : depending on the disorder strength, either L c < 
a or R c > ao is the physically relevant quantity. On 
the other hand, the length L c as defined via fc 4) (0) is 
a convenient parameter quantifying the strength of the 
disorder. 

Assuming L c < ao and short distances L < L c the 
mean-squared displacement behaves as 

(u 2 (L))~e(L/L c ) 3 ; (51) 
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the condition (u 2 (L c )) ~ £ 2 replaces the previous con- 
dition (u 2 (R c )) ~ £ 2 valid in the weak pinning regime. 
Going beyond the collective pinning (or Larkin-) length 
L c the displacement grows as 



with p(r) given by (|24|l . On larger scales R c /a < p(r) < 
A we have found a sub-logarithmic law, 



(u 2 (L))~e(L/L c ) 2 ^ 



(52) 



where £f \ is the roughness exponent for the single vortex. 

Let us now connect the single vortex- with the lattice 
pinning regimes at the scales L ~ clq and R ~ a . In 
the weak pinning case (L c > ao), we match (I51f) with the 
logarithmic growth (JTJ valid for ao < p(r) < R c /a, 

(u 2 (r))~e (a /L c ) 3 ln[p(r)/a ], 

and we recover the relation hi[R c /aao] ~ (L c /ao) 3 , in 
agreement with (f 501) above. The equation for the sub- 
logarithmic growth (|49|l valid for p(r) > R c /a can then 
be casted in the form 



(u 2 (r))~e {{a /L c f\n[p{v)/a 



2a 



For intermediate pinning (£ -C L c < ao) we match (|52[) 
(valid up to L ~ ao) to the sub-logarithmic growth 1|49|) 
and obtain the relation 

(u 2 (r))~e (a /L c ) 2C -(ln[p(r)/a ]) 2ff 

describing the growth of the displacement correlator for 
larger distances p(r) > ao. 

Comparing the two cases of weak and intermediate 
pinning in the sub-logarithmic regime we note a slightly 
different dependence on the disorder parameter L c : in 
the weak pinning situation the disorder dependence is 
given by the factor £~ 6cr , whereas for intermediate pin- 

ning the corresponding factor L c 1,2 involves a different 
exponent. The numerical values of the two exponents 
can be calculated within RG; the individual vortex is 
characterized by an isotropic elasticity and £™ 2 w 0.177e 
to first loop order iSi This value is close to the one ob- 
tained for the anisotropic elasticity of the vortex lattice, 
Cd,2 ~ 0.174e. Therefore, the exponent 2Cj™ 2 ~ 1.062 
is similar to the value 6cr w 1.044. However, note that 
the roughness exponent ^f v 2 is subject to corrections by 
higher-loop terms, whereas the value for a remains un- 
changed. 



V. CONCLUSION 

Summarizing, the long-ranged interaction between vor- 
tices as expressed by the dispersive elastic moduli of the 
flux-line lattice strongly reduces the effect of disorder: 
rather than the usual algebraic growth (w 2 (r)) oc r 2 ^ of 
the mean-squared displacement correlator, one finds a 
mere logarithmic dependence in the perturbative regime 
a < p(r) < Rc/a^. 



»(r)\ ~£ 2 ( ln M r )/ Q °n 



2a 



\ln(Rc/aao) J 

here, the scale function p(r) is determined to the preci- 
sion given in l|20|) with a — a' w 0.32. The exponent 
2(7 w 0.348 obtained within the one-loop RG calculation 
remains unchanged by higher loop corrections. The def- 
inition 



(53) 



(« 2 (r))=C 



2 ln[p(r)/q ] 
In [_R c /aa ] 



for the relevant length scale of the disorder landscape 
emerges naturally within the RG framework. 

The mean-squared displacement correlator for the dis- 
ordered vortex lattice then exhibits two separate regimes 
with a well developed translational order involving either 
a logarithmic or a sub- logarithmic growth of (u 2 (r)), a 
first regime at small length scales originating from the 
long-ranged interaction between vortices, and a second 
regime at large scales which is due to the periodicity of 
the lattice (Bragg Glass regime, cf. Eq. J5J). A number 
of experiments have been devised in order to observe this 
type of logarithmic order, e.g., via Bitter decoration^iiS 
or using small angle neutron scattering (SANS)^*^ (we 
note that the Bitter decoration analysis carried out at 
low fields does not probe the dispersive regime discussed 
here). The observation of a well (i.e., quasi- long range) 
ordered lattice then is often interpreted as evidence for 
Bragg-glass order. The above discussion demonstrates, 
that this need not to be the case, however, as the vor- 
tex lattice exhibits this type of logarithmic order on in- 
termediate scales ao < R < A, ao < L < A 2 /ao as well 
(note the 'cigar' shaped geometry of the dispersive regime 
which can reach a large extension along the longitudinal 
direction, i.e., parallel to the magnetic field). 

A further complication of this type of analysis is found 
in the slow relaxation of the vortex system towards the 
proper glassy order: frozen-in at some higher tempera- 
ture T mo it where pinning is irrelevant, the glassy order of 
the vortex lattice first has to establish itself via proper 
relaxation through creep over the relevant barriers. Dur- 
ing the experimental time t, the vortex system then can 
overcome barriers of size T\n(t/to) typically, with 1 a 
proper attempt frequency. With typical temperatures 
of order 50 K and ln(t/t ) w 20, barriers U(R) w 1000 K 
may be overcome, producing relaxed domains of size (R) 
of a few lattice constants typically; hence reaching glassy 
order over large scales is quite a nontrivial task. 
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